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Abstract 

We  consider  the  broadcasting  problem  in  multi-radio  multi-channel  ad  hoc  networks.  The  objective 
is  to  minimize  the  total  cost  of  the  network-wide  broadcast,  where  the  cost  can  be  of  any  form  that  is 
summable  over  all  the  transmissions  (e.g.,  the  transmission  and  reception  energy,  the  price  for  accessing 
a  specific  channel).  Our  technical  approach  is  based  on  a  simplicial  complex  model  that  allows  us  to 
capture  the  broadcast  nature  of  the  wireless  medium  and  the  heterogeneity  across  radios  and  channels. 
Specifically,  we  show  that  broadcasting  in  multi-radio  multi-channel  ad  hoc  networks  can  be  formulated 
as  a  minimum  spanning  problem  in  simplicial  complexes.  We  establish  the  NP-completeness  of  the 
minimum  spanning  problem  and  propose  two  approximation  algorithms  with  order-optimal  performance 
guarantee.  The  first  approximation  algorithm  converts  the  minimum  spanning  problem  in  simplical 
complexes  to  a  minimum  connected  set  cover  problem.  The  second  algorithm  converts  it  to  a  node- 
weighted  Steiner  tree  problem  under  the  classic  graph  model.  These  two  algorithms  offer  tradeoffs 
between  performance  and  time-complexity.  In  a  broader  context,  this  work  appears  to  be  the  first  that 
studies  the  minimum  spanning  problem  in  simplicial  complexes  and  weighted  minimum  connected  set 
cover  problem. 
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I.  Introduction 

Multi-Radio  Multi-Channel  (MR-MC)  wireless  networking  arises  in  the  context  of  wireless 
mesh  networks,  dynamic  spectrum  access  via  cognitive  radio,  and  next- generation  cellular  net¬ 
works  [1,  2].  By  the  use  of  multiple  channels,  adjacent  transmissions  can  be  carried  over 
non-overlapping  channels  to  avoid  mutual  interference.  Furthermore,  each  node,  equipped  with 
multiple  radios,  is  capable  of  working  in  a  full-duplex  mode  by  tuning  the  transmitting  and 
receiving  radios  to  two  non-overlapping  channels. 

The  increasing  demand  for  high  data  rate  and  the  persistent  reduction  in  radio  costs  have 
greatly  stimulated  research  on  MR-MC  networks.  Considerable  work  has  been  done  on  capacity 
analysis  [1],  channel  and  radio  assignment  [3-8],  and  routing  protocols  [4,  5,  7,  9].  Considerable 
work  has  also  been  done  on  broadcasting  in  Single-Radio  Single-Channel  (SR-SC)  networks  [10- 
13].  In  this  technical  report,  we  consider  the  broadcasting  problem  in  MR-MC  ad  hoc  networks. 

A.  Broadcasting  in  SR-SC  Networks 

Broadcasting  is  a  basic  operation  in  wireless  networks  for  disseminating  a  message  containing, 
for  example,  situation  awareness  data  and  routing  control  information,  to  all  nodes.  For  an  SR- 
SC  network,  a  key  question  for  the  network-wide  broadcast  is  which  set  of  nodes  should  be 
selected  to  transmit  during  the  broadcast  such  that  the  total  cost  is  minimized.  The  cases  where 
the  cost  is  energy  consumption  [10,  11],  the  number  of  transmissions  [12],  or  the  overhead 
in  route  discovery  and  management  [13]  have  been  well  studied.  Recently,  a  unified  solution 
to  the  minimum  cost  broadcasting  problem,  where  the  cost  function  can  have  various  forms, 
has  been  proposed  based  on  the  formulation  of  the  so-called  neighborhood  complex  [14].  In 
contrast  to  their  counterparts  in  wired  networks  which  have  polynomial  solutions,  the  broadcast 
problems  for  minimizing  the  energy  consumption  and  the  number  of  transmissions  are  shown 
to  be  NP-complete  in  [11,  12]. 

The  complexity  of  the  broadcasting  comes  from  the  broadcast  nature  of  the  wireless  medium: 
via  an  omnidirectional  antenna,  a  single  transmission  from  one  node  can  reach  all  the  other  nodes 
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within  the  transmission  range  of  this  node,  but  it  may  cause  interference  to  other  nearby  trans¬ 
missions.  This  “node-centric”  nature  of  the  wireless  broadcasting  problem  along  with  the  mutual 
interference  between  concurrent  transmissions  complicates  the  design  of  efficient  broadcasting 
algorithms. 

B.  Broadcasting  in  MR-MC  Networks 

In  an  MR-MC  ad  hoc  network,  such  as  the  DARPA  Wireless  Network  after  Next  (WNaN)  [15], 
each  node  is  equipped  with  multiple  radios  each  operating  on  a  different  channel.  The  intro¬ 
duction  of  multiple  channels  and  multiple  radios  further  complicates  the  design  of  an  efficient 
broadcasting  scheme.  Since  the  number  of  radios  is  usually  smaller  than  the  number  of  channels, 
the  broadcast  scheme  should  decide  not  only  which  nodes  act  as  relays  but  also  for  those  relay 
nodes,  which  channel(s)  should  be  assigned  to  the  transmitting  radio(s).  Given  the  selection  of 
the  relay  nodes,  two  simple  broadcast  schemes  [5]  are:  (i)  transmitting  multiple  copies  of  the 
message  on  all  channels;  (ii)  transmitting  a  single  copy  of  the  message  on  a  common  channel 
dedicated  to  broadcasting.  Both  schemes  are  inefficient.  For  the  latter  one,  if  the  broadcast  load 
is  high,  the  common  channel  will  be  overwhelmed,  even  while  there  are  plenty  of  other  channels 
free. 

One  subtle  issue  is  the  complication  of  the  wireless  broadcast  advantage.  In  an  MR-MC 
network,  if  the  radios  of  the  neighboring  nodes  are  tuned  to  different  channels,  a  single  trans¬ 
mission  on  one  channel  cannot  reach  all  the  neighboring  nodes  simultaneously.  In  other  words, 
only  the  neighboring  nodes  on  the  same  channel  can  share  the  wireless  broadcast  advantage. 
More  precisely,  the  concept  of  neighborhood  must  be  defined  both  by  radio  range  and  channel. 
Another  subtle  issue  is  channel  heterogeneity.  Channels  may  have  different  bandwidth,  fading 
condition,  and  accessing  cost,  leading  to  different  implications  in  the  total  broadcast  cost. 

Broadcasting  in  MR-MC  networks  is  thus  a  multi-faceted  problem,  involving  channel  assign¬ 
ment,  relay  node  selection,  and  channel  selection  for  the  source  and  relay  nodes.  In  this  technical 
report,  we  focus  on  the  latter  two  issues  by  assuming  a  given  channel-to-radio  assignment.  To 
avoid  the  hidden  channel  problem  [8],  two  nodes  that  are  two-hops  away  from  each  other  are 
assigned  two  distinct  sets  of  channels.  Our  design  objective  is  to  minimize  the  total  broadcast 
cost,  where  the  cost  can  be  of  any  form  that  is  summable  over  all  the  transmissions,  including, 
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for  example,  the  transmission  and  reception  energy1,  the  price  for  accessing  each  channel. 

C.  A  Simplicial  Complex  Model  for  Broadcasting  in  MR-MC  Networks 

Our  technical  approach  is  based  on  a  simplicial  complex  model  of  the  broadcasting  problem  in 
MR-MC  networks.  A  simplicial  complex  is  a  collection  of  nonempty  sets  with  finite  size  that  is 
closed  under  the  subset  operation.  In  other  words,  if  a  set  s  belongs  to  the  collection,  all  subsets 
of  s  also  belongs  to  the  collection.  An  element  of  the  collection  is  called  a  simplex  or  face.  This 
constraint  is  often  satisfied  in  the  network  context.  For  example,  subsets  of  a  broadcast/multicast 
group  are  broadcast/multicast  groups,  subsets  of  a  clique  are  cliques.  A  simple  example  of  graph 
and  simplicial  complex  is  given  in  Fig.  1.  While  the  concept  of  simplicial  complex  has  been 
around  since  the  1920’s,  many  well-solved  fundamental  problems  in  graph  remain  largely  open 
under  this  more  general  model. 


Vo 


Vo 


Vi 


V2 


Vi 


v2 


(a)  Graph 


(b)  Simplicial  Complex 


Fig.  1.  Graph  and  simplicial  complex:  V  =  {vo,vi,V2},  <S(a)  =  {(no,  Vi),  (vo,  v2),  (vi,  u2)},  <$((,)  = 

{(v0,vi,v2),  (v0,vi),  (v0,  v2),  (vi,  v2),  (n0),  (m),  (n2)}. 


We  use  a  simplicial  complex  model  rather  than  a  graph  because  the  simplicial  complex 
more  naturally  captures  the  broadcast  channel,  and  the  distinction  and  disjointness  between 
broadcasting  on  different  channels.  Further,  costs  can  be  attached  to  faces  (simplices)  in  a  way 
not  easily  possible  with  graphs. 

’The  ‘reception  energy’  denotes  the  energy  consumed  by  the  radio  in  reception  mode. 
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Consider  an  example  MR-MC  network.  As  shown  in  Fig.  2,  after  the  channels  are  assigned, 
the  network  is  partitioned  into  cliques  of  nodes.  A  clique  consists  of  the  nodes  which  share  at 
least  one  common  channel,  and  two  cliques  are  spliced  via  nodes  operating  on  multiple  channels 
commonly  shared  by  the  two  cliques.  Within  each  clique,  depending  on  the  cost  function,  the 
transmitter  decides  which  dimension  simplex  (i.e.,  a  subclique  or  the  clique  itself)  in  a  clique 
complex  to  activate.  The  message  for  the  network-wide  broadcast  is  thus  propagated  through  a 
sequence  of  cliques,  possibly  of  different  dimensions.  Note  that  the  unicast  case  corresponds  to 
a  clique  of  dimension  1  (an  edge).  This  example  could  also  apply  to  the  case  where  nodes  may 
have  multiple  radios,  perhaps  of  different  modality  (e.g.,  RF  and  optical);  in  this  case,  there  may 
also  be  a  cost  associated  with  switching  modes. 


Fig.  2.  An  illustration  of  an  MR-MC  network  and  the  constructed  simplicial  complex.  The  parameters  within  the  braces  are  the 
channels  which  each  node  can  access.  In  the  communication  graph  derived  from  the  network,  a  link  exists  between  two  nodes 
if  and  only  if  two  nodes  are  within  each  other’s  transmission  range  and  they  share  at  least  one  common  channel.  Notice  that  a 
clique  in  the  communication  graph  may  not  be  a  clique  in  the  MR-MC  network  (correspondingly,  a  simplex  in  the  simplicial 
complex),  e.g.,  the  three  nodes  of  the  right  empty  triangle. 


In  this  case,  the  network-wide  broadcast  problem  can  be  formulated  as  the  minimum  spanning 
problem  in  simplicial  complexes.  A  clique  in  the  MR-MC  network  is  modeled  as  a  simplex 
in  the  simplicial  complex  (see  Fig.  2),  and  since  a  subset  of  a  clique  is  still  a  clique,  the 
constructed  simplicial  complex  meets  the  requirement  of  being  closed  under  the  subset  operation. 
The  minimum  spanning  problem  in  a  simplicial  complex  is  to  find  a  connected  subset  of 
simplices  that  covers  all  the  vertices  with  the  minimum  total  weight,  i.e.,  the  Minimum  Connected 
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Spanning  Subcomplex  (MCSSub)2.  Then  the  solution  to  the  network-wide  broadcast  problem 
can  be  obtained  by  solving  the  MCSSub  problem. 

D.  Minimum  Spanning  Problem  in  Simplicial  Complexes 

The  minimum  spanning  problem  in  a  graph  is  to  find  a  connected  subgraph  that  covers  all  the 
vertices  with  minimum  total  weight.  The  solution  must  be  a  tree  for  graphs  with  nonnegative 
weights  (hence  called  the  Minimum  Spanning  Tree  (MST)).  There  are  several  polynomial-time 
algorithms  for  MST,  e.g.,  Kruskal’s  Algorithm  and  Prim’s  Algorithm  [16].  The  MST  problem 
has  many  applications  in  network  planning,  broadcasting  in  communication  networks,  touring 
problems,  and  VLSI  design  [17]. 

With  the  addition  of  high  dimensional  simplices,  the  minimum  spanning  problem  in  a  simpli¬ 
cial  complex  is  fundamentally  different  and  much  more  difficult  than  its  counterpart  in  a  graph. 
First,  unlike  the  case  in  a  graph,  the  MCSSub  of  a  simplicial  complex  may  not  be  a  “tree”3.  As 
illustrated  in  Fig.  3,  the  MCSSub  of  the  simplicial  complex  is  the  three  filled  triangles  which 
form  a  cycle.  Second,  while  simple  greedy-type  polynomial-time  algorithms  exist  for  finding  the 
minimum  spanning  tree  in  a  graph,  the  minimum  spanning  problem  in  a  simplicial  complex  is 
NP-complete  as  established  in  this  technical  report  (see  Sec.  III-A). 

We  thus  develop  polynomial-time  approximation  algorithms  for  the  minimum  spanning  prob¬ 
lem  in  simplicial  complexes.  We  propose  two  algorithms:  one  reduces  this  problem  to  a  minimum 
connected  set  cover  problem,  and  the  other  reduces  the  problem  to  a  node-weighted  Steiner 
tree  problem  in  a  graph  derived  from  the  original  simplicial  complex.  We  also  establish  the 
approximation  ratios  of  the  two  algorithms.  Both  are  shown  to  be  order-optimal.  The  time- 
complexity  of  these  two  algorithms  is  also  analyzed,  illustrating  the  tradeoff  between  performance 
and  complexity  offered  by  these  two  algorithms.  In  a  broader  context,  this  work  appears  to  be  the 
first  that  studies  the  minimum  spanning  problem  in  simplicial  complexes  and  weighted  minimum 
connected  set  cover  problem. 

2Strictly  speaking,  a  subcomplex  should  also  be  closed  under  the  subset  operation,  but  without  loss  of  generality,  we  do 
not  include  this  condition  in  the  definition  of  minimum  connected  spanning  subcomplex,  which  is  also  more  relevant  to  the 
broadcasting  problem  at  hand. 

’Although  there  is  no  unified  definition  of  tree  in  simplicial  complexes,  a  couple  of  definitions  can  be  obtained  by  generalizing 
those  equivalent  definitions  of  tree  in  a  graph.  For  example,  simplicial  trees  can  be  defined  based  on  the  universal  existence  of 
leaves  in  any  subgraph,  or  the  uniqueness  of  simplicial  facet  paths. 
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Fig.  3.  A  simplicial  complex  where  its  MCSSub  is  the  three  filled  triangles,  and  is  not  a  “tree”  (the  integers  are  the  weights 
of  the  simplices). 


E.  Related  Work 

Broadcasting  in  MR-MC  networks,  mostly  in  the  context  of  wireless  mesh  networks,  has 
been  studied  for  different  optimization  objectives  [8,  18-23].  A  channel  assignment  algorithm 
that  minimizes  interference  among  nodes  in  a  multicast  routing  tree  is  proposed  in  [8].  Several 
heuristic  algorithms  for  minimizing  the  broadcast  latency  are  presented  in  [18,  19],  where  the 
broadcast  latency  is  defined  as  the  maximum  delay  between  the  transmission  of  a  packet  by  the 
source  node  and  its  eventual  reception  by  all  the  other  nodes.  Efficient  broadcast  algorithms  are 
designed  to  minimize  the  total  number  of  transmissions  in  [20]  or  the  total  number  of  radios 
used  in  the  broadcast4  in  [21].  In  [18-21],  the  NP-hardness  of  the  corresponding  broadcast 
problems  is  established.  The  broadcast  latency  and  the  number  of  transmissions  are  minimized 
simultaneously  in  [22].  A  widest  spanning  tree  is  constructed  in  [23]  to  maximize  the  bottleneck 
link  bandwidth  of  the  tree  such  that  the  content  update  time  is  minimized. 

4Notice  that  the  total  number  of  radios  counts  not  only  the  transmitting  radios  but  also  the  receiving  radios,  and  thus  it  is 
strictly  greater  than  the  total  number  of  transmissions. 
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Different  from  the  previous  work,  the  optimization  objective  in  our  work  can  be  any  cost 
function  which  is  summable  over  all  the  transmissions,  thus  taking  into  account  channel  het¬ 
erogeneity  (e.g.,  transmissions  on  different  channels  may  consume  different  amounts  of  energy, 
due  to  different  bandwidths  or  different  propagation  characteristics  or  some  other  factor).  We 
point  out  that  neither  minimizing  the  total  number  of  transmissions  nor  minimizing  the  total 
number  of  radios  used  in  the  broadcast  is,  in  general,  equivalent  to  minimizing  the  total  energy 
consumption.  The  reception  energy  is  ignored  if  the  total  number  of  transmissions  is  minimized, 
while  the  transmission  energy  and  the  reception  energy  are  equated  if  the  total  number  of  radios 
is  minimized.  More  importantly,  channel  heterogeneity  is  not  addressed  if  these  two  objectives 
are  optimized. 

Furthermore,  to  our  best  knowledge,  our  work  is  the  first  to  adopt  simplicial  complexes  to 
model  and  solve  the  broadcast  problem  in  wireless  ad  hoc  networks.  An  attempt  to  formulating 
the  broadcast  problem  in  wireless  ad  hoc  networks  into  problems  in  simplicial  complexes  can 
be  found  in  [14].  The  simplicial  complex  is  an  important  topic  in  algebraic  topology.  A  detailed 
analysis  based  on  homology  and  cohomology  can  be  found  in  [24,  25].  Homology  theory  and 
techniques  have  been  adopted  as  a  computation  tool  for  detecting  and  locating  the  coverage 
holes  in  sensor  networks  [26-28],  and  for  deducing  the  topological  information  of  the  physical 
space  from  the  encounter  data  of  mobile  networks  [29]. 

Simplicial  complexes  can  also  be  considered  as  special  hypergraphs  with  the  added  constraint 
on  the  closedness  under  the  subset  operation.  Here  we  adopt  the  simplicial  complex  model 
instead  of  the  more  general  hypergraph  model,  because  broadcasting  does  exhibit  subset  closure 
(, i.e .,  one  node  can  choose  to  broadcast  to  any  subset  of  other  nodes  within  the  same  clique). 

II.  Basic  Concepts  in  Simplicial  Complexes 

In  this  section,  we  introduce  several  basic  concepts  in  simplicial  complexes  [24,  25]. 

An  (abstract)  simplicial  complex  is  a  collection  A  of  nonempty  sets  with  finite  size  such  that 
if  A  e  A,  then  V  B  C  A,  B  e  A,  i.e.,  A  is  closed  under  the  operation  of  taking  subsets.  The 
element  A  of  A  is  called  a  simplex  of  A;  its  dimension  (denoted  by  dim  A)  is  one  less  than  the 
number  of  its  elements.  Each  nonempty  subset  of  A  is  called  a  face  of  A.  The  dimension  of  A 
is  the  maximum  dimension  over  all  its  simplices,  or  is  infinite  if  the  maximum  does  not  exist. 
The  vertex  set  V  of  A  is  the  union  of  the  one-point  elements  of  A.  Fig.  4  shows  an  example 
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of  a  2-dimensional  simplicial  complex.  A  subcollection  of  A  that  is  itself  a  simplicial  complex 
is  called  a  subcomplex  of  A.  A  subcomplex  of  A  is  the  p-skeleton  of  A,  denoted  by  A if  it 
is  the  collection  of  all  simplices  of  A  with  dimension  no  larger  than  p.  Thus,  the  1-skeleton  is 
the  underlying  graph  of  A. 


Fig.  4.  A  simplicial  complex  A  with  6  vertices  (0-dimensional  simplices:  {-t!o},{'t>i},{t;2},{i’3},{'t'4},{v5}),  5  edges  (1- 
dimensional  simplices:  {i>i,  V2},{v2,  V3 },{«3 ,  V4},{r>3,  vs},{v4,  V5}),  and  1  filled  triangle  (2-dimensional  simplex:  {113,1:4,115}). 


Any  abstract  simplicial  complex  has  a  geometric  realization.  There  is  a  one-to-one  corre¬ 
spondence  between  each  element  A  =  {u0, vn}  of  A  and  a  simplex  spanned  by  the  vertices 
Vo,...,vn  in  Euclidean  space.  It  can  be  shown  that  every  d-dimensional  simplicial  complex  can 
be  realized  in  M2d+1.  Take  graphs  as  a  special  example,  which  are  1-dimensional  simplicial 
complexes:  non-planar  graphs  cannot  be  embedded  in  M2,  but  all  the  graphs  can  be  realized  in 
M3. 

A  facet  of  a  simplicial  complex  A  is  a  maximal  face  of  A,  i.e.,  it  is  not  a  subset  of  any  other 
face.  A  simplicial  complex  is  connected  if  its  1-skeleton  (i.e.,  the  underlying  graph)  is  connected 
in  the  graph  sense. 

A  weighted  simplicial  complex  (WSC)  A  is  a  triple  (V,S,  w)5,  where  V  is  the  set  of  vertices, 
S  the  set  of  faces  of  A,  and  w  :  S  — *  (M+  U  {0}}  a  nonnegative  weight  function  defined  for 
each  face  in  S  with  w(v)  =  0  for  all  v  G  V.  We  define  the  facet-only  weight  WF(  A)  of  a  WSC 


5(S,w)  suffices  to  denote  the  WSC  since  V  C  S,  but  we  use  the  redundant  (V,S,w)  for  convenience. 
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A  as 

WF{  A)  =  ^  w(Ft). 

FiC {facet  of  A} 

III.  Minimum  Connected  Spanning  Subcomplex 

In  this  section,  we  show  that  the  MCSSub  problem  is  NP-complete,  and  we  propose  two 
approximation  algorithms  based  on  connected  set  cover  and  node-weighted  Steiner  tree.  We  also 
establish  the  approximation  ratios  of  the  two  algorithms  and  analyze  their  time  complexity. 

A.  NP-Completeness 

The  decision  version  (D-MCSSub)  of  the  MCSSub  problem  is  stated  as  follows:  let  C(A) 
denote  the  vertex  set  of  a  WSC  A  and  WF{ A)  the  facet-only  weight  of  A.  Given  a  WSC 
A  =  (V.  S,  w)  and  K  >  0,  is  there  a  connected  subcomplex  Asub  of  A  such  that  V ( Asub )  =  V 
and  WF{ Asub)  <  K1  Then  we  have  the  following  theorem. 

Theorem  1:  The  D-MCSSub  problem  is  NP-complete. 

To  prove  the  NP-completeness,  we  reduce  the  classic  NP-complete  problem  -  the  unweighted 
set  cover  problem  to  the  MCSSub  problem. 

Proof:  To  check  a  solution  to  the  D-MCSSub  problem,  we  only  need  to  verify  the  following 
points:  (i)  compute  the  facet-only  weight  of  the  solution  and  compare  the  weight  with  K ;  (ii) 
check  whether  the  underlying  graph  of  the  solution  is  connected  or  not;  (iii)  check  whether  all 
the  vertices  of  the  original  simplicial  complex  are  covered  by  the  solution.  Since  all  these  can 
be  done  within  polynomial  time,  the  D-MCSSub  problem  is  NP. 

Given  is  an  unweighted  Set  Cover  instance  I,  i.e.,  a  universe  of  elements  U  and  a  family 
of  subsets  F  of  U.  For  each  element  u  £  U,  we  introduce  a  corresponding  vertex  in  MCSSub 
instance  V .  We  introduce  one  additional  vertex  d.  For  each  set  /  e  F,  we  introduce  a  corre¬ 
sponding  face  f  —  f  U  {d},  of  weight  1.  Being  a  simplcial  complex,  all  subsets  of  /'  are  also 
introduced,  all  also  of  weight  1.  In  addition  to  covering  all  vertices,  a  solution  to  I'  must  be 
connected. 

Given  any  solution  SOL  to  /  of  cost  c,  we  can  construct  a  solution  to  P  also  of  cost  c.  For 
every  set  /  in  SOL  not  including  d,  replace  /  with  /  U  {d}. 

On  the  other  hand,  given  any  solution  SOL'  to  /'  of  cost  c,  a  solution  to  /  of  the  same  cost 
can  be  constructed  as  follows:  for  any  face  /  6  SOL'  such  that  /  does  not  appear  as  a  set  in 
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F,  replace  /  with  any  superset  of  f  —  {d}  appearing  in  F.  Note  that  there  must  exist  at  least 
one  such  superset.  ■ 

This  proof  is  for  general  D-MCSSub  problems.  It  can  be  shown  that  even  if  the  weight  function 
of  the  WSC  is  monotone  or  strictly  monotone6,  the  D-MCSSub  problem  is  NP-hard.  But  it  is 
still  possible  that  the  D-MCSSub  problem  under  some  special  structured  weight  function  is  P. 

In  the  following,  we  present  two  approximation  algorithms  for  the  MCSSub  problem  both 
with  performance  guarantee  O(lnn),  where  n  is  the  number  of  vertices  in  the  WSC.  Since  the 
best  possible  approximation  ratio  for  the  set  cover  problem  is  Inn  [30],  these  two  algorithms 
are  order-optimal. 

B.  Algorithm  Based  on  Connected  Set  Cover 

Let  A  be  a  set  with  finite  number  of  elements,  and  B  =  {f?*  C  A  :  i  —  1, ...,  n}  a  collection 
of  subsets  of  A  where  each  Bi  is  associated  with  a  weight  w(Bj)  >  0.  Let  G  be  a  connected 
graph  with  the  vertex  set  B.  A  connected  set  cover  (CSC)  Sc  with  respect  to  (A,  B,  w,  G )  is  a 
set  cover  of  A  such  that  Sc  induces  a  connected  subgraph  of  G.  The  minimum  connected  set 
cover  (MCSC)  problem  is  to  find  the  CSC  with  the  minimum  weight,  where  the  weight  of  a 
CSC  Sc  is  defined  as 

IU(SC)  =  ^2 

BiCSc 

From  a  WSC  A  =  (V,S,w),  we  derive  an  auxiliary  undirected  graph  Ga  in  the  following 
way:  let  S  \  V  be  the  vertex  set  of  Ga,  and  connect  two  vertices  (non-vertex  faces  in  A)  Si 
and  S2  if  and  only  if  S\  ft  S2  f  0.  Then  we  have  the  following  theorem  on  the  relation  between 
the  MCSSub  problem  and  the  MCSC  problem. 

Theorem  2:  Let  A*  be  the  MCSSub  of  a  WSC  A  =  (V.  S ,  w)  and  Sf  the  MCSC  of  ( V ,  S  \ 
V,w,G a)-  Then  we  have 

wf(A*)  =  w(Sc). 

Proof:  The  proof  is  based  on  the  following  lemma. 

Lemma  1:  Let  Sf  be  the  MCSC  of  (V,  S\V,w,  Ga)-  For  any  face  S  G  Sf  with  w(S)  >  0, 
we  have  that  there  does  not  exist  a  face  S'  E  S*,  such  that  S  C  S'. 

6We  say  that  the  weight  function  satisfies  the  monotone  property  if  for  any  two  faces  Si  C  S2,  w(Si)  <  w(S2),  i.e.,  the 
weight  is  monotone  non-decreasing  with  respect  to  the  dimension  of  the  face. 
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Proof  of  Lemma  1:  Suppose  that  for  some  face  S  G  Sf  with  w(S)  >  0,  3  S'  E  Sf  such  that 
S  C  S'.  Let  S'c  =  Sf\  s.  Obviously,  S'c  is  a  set  cover,  and  w(S'c)  =  w(Sf)  —  w(S )  <  w(Sf). 
On  the  other  hand,  since  S  D  S"  0  implies  S'  Cl  S”  f  0  for  any  face  S"  G  Sf,  it  follows  from 
the  connection  rule  of  the  auxiliary  graph  G a  that  any  path  via  S  has  an  alternative  path  via  S'. 
Thus,  S’c  is  a  CSC,  leading  to  a  contradiction.  ■ 

Given  the  MCSC  Sf  of  {V,  S\V,w,  Ga),  we  can  obtain  a  connected  spanning  subcomplex 
by  mapping  each  element  of  Sf  to  a  face  in  A.  Since  the  facet-only  weight  wf(^c)  °f 
only  counts  facets  in  Ajf,,  it  follows  that  iup(A *c)  <  w(Sf).  Based  on  Lemma  1,  we  have  that 
every  element  of  S*c  with  positive  weight  is  a  facet  in  Af,  and  thus 

wf( A*)  <  wf(A*c )  =  w(Sf), 

where  A*  is  the  MCSSub  of  A. 

On  the  other  hand,  the  facets  of  A*  leads  to  an  CSC  S^,  and  w(S^f)  =  wF( A*).  It  implies 
that 

w{S*c)  <  w(S*A)  =  wf(A*). 

Thus,  wf(A*)  =  w(Sf).  M 

1)  Algorithm:  Based  on  Theorem  2,  we  can  reduce  the  MCSSub  problem  of  a  WSC  A  = 
(V.S,w)  to  the  MCSC  problem  (L,  S  \  V,w,GA).  We  obtain  the  following  Set  Cover  based 
Algorithm  (SCA)  for  the  MCSSub  problem. 


Algorithm  1:  SCA  for  MCSSub: 

INPUT:  A  WSC  A  =  (V,S,w). 

OUTPUT:  An  approximate  MCSSub  Ac  of  A. 

1.  Derive  the  auxiliary  graph  G A. 

2.  Find  an  approximate  MCSC  Sc  of  ( V.  S  \  V,w,GA )  by  using  the  greedy  algorithm  for 
MCSC  (Algorithm  2). 

3.  Transform  Sc  to  a  connected  spanning  subcomplex  Ac  by  mapping  each  element  of  Sc  to 


a  face  in  A. 
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Zhang  et  al.  propose  a  greedy  approximation  algorithm  for  the  unweighted  MCSC  prob¬ 
lem  [3 1]7,  i.e.,  w(Bj)  =  1  for  all  i.  By  generalizing  their  greedy  approach,  we  develop  a  greedy 
algorithm  for  the  weighted  MCSC  problem. 

Before  stating  the  algorithm,  we  introduce  the  following  notations  and  definitions.  For  two 
sets  3i ,  Sr2  G  S,  let  distf;(Sj ,  S2)  be  the  length  of  the  shortest  path  between  Sj  and  S2  in  an 
auxiliary  graph  G,  where  the  length  of  a  path  is  given  by  the  number  of  edges;  Si  and  S2  are 
said  to  be  graph-adjacent  if  they  are  connected  via  an  edge  in  G  (i.e.,  distG(S'i,  S2)  =  1),  and 
they  are  said  to  be  cover-adjacent  if  Si  D  S2  ^  0.  Notice  that  in  a  general  MCSC  problem,  there 
is  no  connection  between  these  two  types  of  adjacency.  The  c over- diameter  Dc(G)  is  defined 
as  the  maximum  distance  between  any  two  cover- adjacent  sets,  i.e., 

DC(G )  =  max{distc(Sj .  S2)  \  Si,  S2  G  S  and  Si  n  S2  0}. 

For  the  MCSC  problem  derived  from  the  MCSSub  problem  of  a  WSC  A,  we  have  that  Dc(GA )  = 

1. 

At  each  step  of  the  algorithm,  let  1Z  denote  the  collection  of  the  subsets  (faces  of  A)  that 
have  been  selected,  and  U  the  vertex  subset  of  A  that  has  been  covered.  Given  1Z  0  and  a  set 
S  E  S\7Z,  an  1Z  — >  S  path  is  a  path  {S'0,  Si, ...,  Sk}  in  G  such  that  (i)  S0  G  1Z\  (ii)  Sk  =  S; 
(iii)  Si, ...,  Sk  G  S  \  JZ.  We  define  the  weight  ratio  r(Ps )  of  Ps  as 

(P  )  _  w(S(ps)  \fc)  _  ^ses(Ps)\nw(S) 

1  |Vat(Ps)|  |V)v(Ps)|  ’  1  ; 

where  S(Ps)\1Z  is  the  subsets  (faces  in  S)  of  Ps  that  are  not  in  7Z,  and  |Cv(Ps)|  is  the  number 
of  vertices  of  A  that  are  covered  by  Ps  but  not  covered  by  7Z. 


Algorithm  2:  A  Greedy  Algorithm  for  MCSC. 

INPUT:  (U  S\V,w,  GA) 

OUTPUT:  A  CSC  1Z. 

1.  Choose  So  G  S  \  V  such  that  the  weight  ratio  r(S0 )  defined  in  (1)  is  the  minimum,  and 
let  1Z  =  {S0}  and  U  =  S0. 

2.  WHILE  V  \  U  ±  0  DO 

7The  original  algorithm  in  [31]  has  a  flaw  and  the  established  approximation  ratio  is  incorrect.  In  [32],  the  flaw  is  corrected 
and  a  stronger  result  on  the  approximation  ratio  is  shown. 
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2.1.  For  each  S  e  S\(VUlZ)  which  is  cover- adjacent  or  graph- adjacent  with  a  set  in  72, 
find  a  shortest8  72  — *  S  path  Ps. 

2.2.  Select  Ps  with  the  minimum  weight  ratio  r(Ps)  defined  in  (1),  and  let  72  =  72  U  Ps 
(add  all  the  subsets  of  Ps  to  72)  and  17=17  U  \fi(Ps). 

END  WHILE 
3.  RETURN  72. 


2)  Approximation  Ratio:  The  approximation  ratio  of  SCA  is  determined  by  Step  2,  i.e.,  the 
approximation  ratio  of  the  greedy  algorithm  for  the  MCSC  problem.  First,  we  show  the  following 
lemma. 

Lemma  2:  Given  a  weighted  MCSC  problem  (V,  S\V,w,  G )  with  Dc(G)  =  1,  let 


H/r 


maxlufiS')) 

Ses  1  v 

nun-furs')} 

ScS  1  J 


(2) 


Then  the  approximation  ratio  of  the  greedy  algorithm  for  MCSC  is  at  most  Rw  +  Ilf:  —  1), 
where  7  =  max{ | Sj  |  S  G  S  \  V}  is  the  maximum  size  of  the  subsets  in  S  and  H(-)  is  the 
harmonic  function. 

Proof  of  Lemma  2:  The  proof  is  based  on  the  classic  charge  argument.  Let  S*  be  an 
optimal  solution  to  the  weighted  set  cover  problem  (V,«S  \  V,w),  and  72  the  solution  returned 
by  the  greedy  algorithm  for  the  weighted  MCSC  problem  (V,  S\V,w,  G )  with  DC(G )  =  1.  Let 
w{S*)  and  w{72)  denote  the  total  weight  of  the  subsets  included  in  S*  and  72,  respectively.  In 
the  following,  we  will  show  that 


^-<^  +  ^(7-1).  (3) 

w(o*) 

Let  Sf  be  an  optimal  solution  to  the  weighted  MCSC  problem  (V,  S\V,  w,  G)  with  DC(G )  = 
1.  Since  w(S*)  <  w(Sf),  Lemma  2  follows  immediately. 

To  prove  (3),  we  apply  the  classic  charge  argument.  Each  time  a  subset  So  (at  step  1)  or 
a  shortest  72  — >  S  path  If  (at  step  2)  is  selected  to  be  added  to  72,  we  charge  each  of  the 
newly  covered  elements  (at  step  1)  or  r(Pg )  defined  in  (1)  (at  step  2).  Notice  that  when 


'Notice  that  the  shortest  path  is  defined  in  terms  of  the  number  of  edges,  not  the  total  weight  of  all  vertices  along  the  path. 
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DC(G )  =  1,  the  shortest  1Z  — >■  S  path  PJ  is  only  a  single  edge  connecting  some  subset  in  1Z 
and  S,  and 


r(P*s)  = 


w(S(P§\K))  _  w(S) 


\VN(P*)\  \s\u\- 

During  the  whole  procedure,  each  element  of  V  is  charged  exactly  once.  Assume  that  step  2 
is  completed  in  K  —  1  iterations.  Let  P^t  be  the  shortest  JZ  — >  S  path  selected  by  the  algorithm 
at  iteration  i.  Let  C'(v)  denote  the  charge  of  an  element  v  in  V.  Then  we  have  that 

E  E  cm 

*=0  vCVN{P*Si ) 


Ec( 

vev 


v  = 


K- 1 

E  E 

*=o  veSi\u 

K- 1 


w(Si) 
\Si\U\ 

^w(Si)  =  w(K), 

i= 0 


(4) 


where  E*0  =  { ,S'0 } . 

Suppose  that  S*  =  {S'J', ...,  S^}  is  a  minimum  weighted  set  cover  for  {V,  S  \  V,  w}.  Since  an 
element  of  V  may  be  contained  in  more  than  one  subset  of  S*,  it  follows  that 


N 


E  cw  <  E  E  cw-  <5) 

vev  1=1  ves* 

Next  we  will  show  an  inequality  which  bounds  from  above  the  total  charge  of  a  subset  in  5*, 
i.e.,  for  any  B*  e  S*, 

C(v)  <  lR™  +  H{\S*\  -  l)]w(S*).  (6) 

v€lB* 

Let  rii  (i  —  0, 1, ...,  K)  be  the  number  of  elements  of  S*  that  have  not  been  covered  by  S  after 
iteration  i  —  1,  where  step  1  is  considered  as  iteration  0.  Let  { i \ .  ...,4}  denote  the  subsequence 
of  {*  =  0, 1, ...,  K  —  1}  such  that  n*  —  ni+ 1  >  0.  For  each  element  a  covered  at  iteration  4,  if 
4  =  0,  based  on  the  greedy  rule  at  step  1,  we  have  that 

w(S*) 


C(v)  =  r(PL )  < 


n, 


(7) 


Otherwise, 


c(v)  =  r(P*  )  =  w(sh)  <  w(S*)R*, 

Sii  \  C/|  nil-n(il+ 1)' 


(8) 
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The  inequality  in  (8)  is  due  to  the  fact  that  Sn  covers  at  least  nn  —  n^1+ 1)  elements  of  V,  i.e., 
|*5u  \U\>  nh  -  n(i1+ 1).  Summing  up  (7)  and  (8), 

w(S*)Rw 


C(v)< 


n. 


h  ™(ii+l) 


(9) 


Consider  two  cases: 

(i)  If  all  the  elements  of  S*  have  been  covered  after  iteration  i\,  i.e.,  n^1+ 1)  =  0,  then 

y  c(v)  <  y  =  w^s*)rw.  ao) 

n° 


(ii)  If  not  all  the  elements  of  S*  have  been  covered  by  1Z  after  iteration  i\,  S*  becomes  cover- 
adjacent  with  1Z  and  thus  a  candidate  for  being  selected  at  the  following  iterations.  At  each 
iteration,  for  each  element  v  G  S*  covered  at  iteration  i3  (j  =  2, k),  the  greedy  rule  at 
step  2  still  yields 


C(v)  =  r(PJJ<r(P| .) 

x3 

w(S*)  _  w(S*) 

WW\~  ' 


(ii) 


It  follows  from  (9,11)  that 


yC'(u)  <  {nh-n{il+ 1))— — 


ves* 


-n{ij+1))  — 
3=  2  lj 

k 

rii.  —  rii,.,,, 

h  lO+i) 


3= 2 


nU 


(12) 


Here  we  have  used  the  fact  that  n^ij+ 1)  =  n,;(j+1) .  It  is  because  between  iteration  i3  and 
iteration  i(j+\),  no  elements  of  S*  are  covered. 

For  the  summation  term  in  (12),  we  have  the  following  inequality: 

1 


nh  niu+i)  <  1  j 

rk,  ~  z--'  n, 


£ 

3= 2 


+ 


ni0.+D  +  1 


J=2 

=  H(ni2)<H(\S*\-l). 


(13) 


The  last  inequality  is  due  to  the  fact  that  ni2  <  —  1  =  \S*\  —  1. 
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Eqn.  (6)  is  a  direct  consequence  of  (10),  (12),  and  (13).  Thus,  using  (4-6), 


N 


w(R)  =  £C(„)<££C(, 


vev 

N 


i=  i  vest 


2—1 


<  [Rw  +  H('y-l)]w(S*). 


Then,  as  a  direct  consequence  of  Lemma  2,  we  have  the  following  theorem  on  the  approxi¬ 
mation  ratio9  of  the  greedy  algorithm  for  the  MCSC  problem  with  Dq(G)  =  1. 

Theorem  3:  Let  A*  be  the  MCSSub  of  a  WSC  A  =  ( V, ,  S ,  w)  and  Ac  be  the  solution  returned 
by  Algorithm  1.  Let  Rw  be  defined  as  in  (2).  Then  we  have 

wf(A  c) 


wf(A* 


<  Rw  +  f/(dimA), 


where  dimA  is  the  dimension  of  A  and  H(-)  is  the  harmonic  function. 

Lrom  Theorem  3,  we  see  that  the  approximation  ratio  depends  on  the  ratio  Rw  of  the  maximum 
weight  to  the  minimum  weight.  It  is  shown  in  the  following  theorem  that  if  Rw  is  unbounded, 
then  the  scaling  order  of  the  approximation  ratio  can  be  as  bad  as  linear  with  respect  to  the 
number  of  vertices  in  the  simplicial  complex. 

Theorem  4:  Let  n  be  the  number  of  the  vertices  in  a  WSC  A  =  (V,S,w),  and  R.w  defined 
as  in  (2).  If  Rw  is  unbounded,  then  the  approximation  ratio  of  Algorithm  1  for  the  MCSSub 
problem  of  A  is  f2(n). 

Proof:  Consider  a  specific  example:  A  is  a  (n  —  1) -dimensional  simplex  with  the  vertex  set 
V  =  {ui,  ...,vn},  and  all  the  weights  of  the  faces  are  infinite  except  for  the  following  five  faces: 


1 

2’ 

1 

"  2’ 


tc(S'i) 

=  w 

({ui,..,ut})  =  i 

w(S2) 

=  w 

({fi, -.•'s, »(,«)}) 

w(S3) 

=  w 

({”(5+1). -.•’(t+i)}) 

w(Sf) 

=  w 

w(S5) 

=  w 

9The  approximation  ratio  of  the  greedy  algorithm  for  general  weighted  MCSC  problem  is  still  an  open  problem. 
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For  ease  of  presentation,  we  have  assumed  that  n  is  a  multiple  of  4.  By  applying  Algorithm  1, 
we  reduce  the  MCSSub  problem  for  A  to  the  MCSC  problem  (V,  S  \  V,w,Ga)-  Due  to  the 
weight  assignment,  it  suffices  to  only  consider  the  subgraph  of  Ga  induced  by  the  above  five 
faces,  as  shown  in  Fig  5. 


Si 


.,V( 


K1) 


Fig.  5.  The  subgraph  of  G a  induced  by  the  five  faces  Si,  S2,  S3,  S4,  and  S5  with  finite  weights. 


The  optimal  solution  A*  to  the  MCSSub  problem  is  given  by 

A*  =  {S  e  S  |  S  C  S2  or  S3  or  S5}, 

and 

wF(  A*)  =  w(S2)  +  iv(S3)  +  w(S5)  =  2. 

On  the  other  hand,  the  solution  Ac  returned  by  Algorithm  1  is  given  by 

Aa  =  {S  e  S  |  S  C  S3  or  ^4}, 


and 


77 

wf(Ac)  =  w(Si)  +  w(S4)  =  -  +  - . 

Z  o 
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Specifically,  Si  is  firstly  selected,  and  then  S4.  Thus, 

Wf(Ac)  =  n_  1  = 

Wp(A*)  16  4 

It  follows  that  the  approximation  ratio  of  Algorithm  1  is  Q(n).  ■ 

From  Theorem  4,  we  see  that  Algorithm  1  is  not  suitable  for  the  MCSSub  problem  of  a  WSC 
A  if  its  weight  function  has  a  relatively  wide  range.  As  shown  next  in  Sec.  III-C,  the  other 
approximation  algorithm  based  on  the  Steiner  tree  does  not  have  this  issue:  its  approximation 
ratio  does  not  depend  on  the  range  of  the  weight  function. 


C.  Algorithm  Based  on  Steiner  Tree 

From  a  WSC  A  =  (V,  S,w),  we  derive  an  undirected  graph  HA  with  the  vertex  set  S:  for 
each  face  S  e  S\V  (i.e.,  the  faces  that  are  not  the  vertices  of  A),  we  replace  it  by  a  vertex  vs 
in  HA  and  connect  vs  to  all  the  vertices  of  S.  The  weight  w(vs)  assigned  to  the  vertex  vs  is 
the  weight  w(S)  of  the  face  S.  Notice  that  the  weight  of  vertices  in  HA  corresponding  to  the 
vertices  in  A  (i.e.,  V)  is  zero.  Fig.  6  shows  an  example  of  the  derivation  of  the  graph  from  a 
2-simplex.  We  have  the  following  theorem  on  the  relation  between  the  MCSSub  of  A  and  the 
Steiner  tree  of  HA  that  spans  the  vertex  set  V  of  A  and  the  minimum  connected  dominating 
set10  of  Ha. 

Theorem  5:  Let  A*  denote  the  MCSSub  of  a  WSC  A  =  ( V. ,  S,  w ),  T*  the  Steiner  tree  of  HA 
that  spans  the  vertex  set  V  of  A,  and  D*c  the  minimum  connected  dominating  set  of  HA.  Then 
we  have  that 

wF(  A*)  =  w(T*)  =  w(D*c). 

Proof:  First  we  show  that  wF(A*)  =  w(T*).  Since  every  connected  spanning  subcomplex 
A'  of  A  corresponds  to  a  connected  subgraph  of  HA  which  only  contains  the  vertices  of  A  and 
the  vertices  representing  the  facets  of  A',  it  follows  that  w(T*)  <  vjf(A*).  On  the  other  hand, 
since  by  contradiction,  there  is  a  one-to-one  mapping  between  the  vertices  of  the  Steiner  tree  of 

10 A  dominating  set  of  a  graph  is  a  subset  of  vertices  such  that  every  vertex  of  the  graph  is  either  in  the  subset  or  a  neighbor 
of  some  vertex  in  the  subset,  and  a  connected  dominating  set  (CDS)  is  a  dominating  set  where  the  subgraph  induced  by  the 
vertices  in  the  dominating  set  is  connected.  The  CDS  problem  asks  for  a  CDS  with  the  minimum  total  weight,  and  it  is  shown 
to  be  a  special  case  of  the  MCSC  problem  [32]. 
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A 


Ha 


Fig.  6.  The  derived  graph  of  a  2-simplex  (squares  in  Ha  represent  the  faces  that  are  not  vertices  of  A). 


Ha  and  the  vertices  plus  the  facets  of  a  connected  spanning  subcomplex  of  A,  it  follows  that 

wF( A*)  <  w(T*). 

Next  we  show  that  w(T*)  =  w(D*c).  Notice  that  the  vertex  set  V  of  A  is  a  dominating  set 
of  Ha-  Since  the  Steiner  tree  T*  of  Ha  spans  the  vertex  set  V,  T*  is  a  CDS  of  Ha-  Thus, 
u:(D*c)  <  w(T*).  On  the  other  hand,  given  the  minimum  CDS  D*c  of  Ha,  since  each  vertex  v 
in  the  vertex  set  V  is  either  in  D*c  or  a  neighbor  of  some  face  in  D*c  and  the  weights  of  the 
vertices  in  V  are  all  zero,  the  combination  of  V  and  I)*c  yields  a  connected  subgraph  of  Ha 
that  spans  V  with  the  same  weight  as  D*c.  Thus,  wiT*)  <  w(D*c).  ■ 

Based  on  Theorem  5,  we  propose  the  following  Steiner  Tree  based  Algorithm  (STA)  for  the 
MCSSub  problem. 


Algorithm  3:  STA  for  MCSSub: 

INPUT:  A  WSC  A  =  (V,S,w). 

OUTPUT:  An  approximate  MCSSub  Ac  of  A. 

1.  Derive  the  graph  Ha  from  A. 

2.  Obtain  an  approximate  Steiner  tree  T  of  Ha  by  using  the  algorithms  given  in  [33,  34]. 

3.  Transform  T  to  a  connected  spanning  subcomplex  A c  of  A  by  mapping  each  element  of 


T  to  a  face  of  A. 
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Since  approximation  only  occurs  in  Step  2,  the  approximation  ratio  of  STA  is  equal  to  that 
of  the  algorithm  for  the  node-weighted  Steiner  tree  problem.  The  best  approximation  ratio  is 
known  to  be  (1.35  +  e)  In  n  for  any  constant  e  >  0,  where  n  is  the  number  of  vertices  of  A  and 
is  also  the  number  of  terminals  in  the  Steiner  tree  of  HA  [34].  Here  we  do  not  try  to  find  the 
CDS  D*c  of  II a.  at  step  2,  because  the  best  known  approximation  ratio  for  the  CDS  problem  is 
(1.35  +  e)  lnn(jfA)  [34,  35].  Since  ri(HA)  n,  the  latter  approximation  ratio  is  much  worse  than 
the  former  one. 

D.  Time  Complexity  Analysis 

Here  we  analyze  the  time  complexity  of  SCA  and  STA  for  the  MCSSub  problem.  Given  a 
WSC  A  =  (V,S,w),  let  n  =  |H|  denote  the  number  of  vertices  in  A,  m  —  |«S  \  V\  the  number 
of  non-vertex  faces  in  A,  and  d  the  dimension  of  A.  Recall  that  the  existence  of  edges  in  the 
auxiliary  graph  GA  for  SCA  and  the  derived  graph  HA  for  STA  depends  entirely  on  whether 
the  two  non-vertex  faces  overlap  and  whether  the  vertex  is  contained  in  the  non-vertex  face, 
respectively.  It  implies  that  all  the  information  of  these  two  graphs  can  be  easily  retrieved  from 
the  WSC  A.  Thus,  Step  1  in  both  algorithms  can  be  skipped  in  the  implementation,  and  the 
time  complexity  of  both  algorithms  is  determined  by  their  Step  2. 

Step  2  of  SCA  is  to  apply  the  greedy  algorithm  to  the  MCSC  problem  (V,S  \  V,  w,  GA).  It 
takes  0(m)  time  to  complete  Step  1  of  the  greedy  algorithm.  Since  at  least  one  vertex  becomes 
covered  at  each  iteration  of  Step  2  of  the  greedy  algorithm,  there  are  at  most  n  —  1  iterations. 
At  each  iteration,  the  weight  ratios  of  at  most  m  faces  are  computed,  and  due  to  the  fact  that 
cover- adjacent  faces  are  graph-adjacent,  the  weight  ratio  of  each  face  is  done  in  constant  time. 
Thus,  the  running  time  of  SCA  is  0(m  +  nm)  =  0(nm). 

Since  the  derived  graph  HA  has  n  +  m  vertices  and  0(dm )  edges  and  the  Steiner  tree  has  n 
terminals  to  cover,  it  follows  from  [36]  that  the  running  time  of  Step  2  of  STA  is  0(dnm 2  + 
nm2  log  m).  From  the  above,  we  see  that  the  time  complexity  of  STA  is  significantly  higher  than 
that  of  SCA.  This  is  mostly  because  the  approximation  algorithm  for  the  Steiner  tree  requires 
the  computation  of  the  shortest  paths  between  all  vertex  pairs. 

We  point  out  that  while  the  Steiner  tree  based  algorithm  has  a  higher  complexity,  it  can  offer 
better  performance  in  a  WSC  with  a  large  weight  range.  In  a  simulation  example  of  random 
simple  complexes,  we  consider  a  case  where  each  face  weight  takes  only  two  values  wmin  and 
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wmax  with  equal  probability.  With  wmin  =  1,  tv-max  =  10000,  and  1000  Monte  Carlo  runs  for 
a  200-vertex  random  simplicial  complex11  [38],  we  find  that  the  total  weight  of  the  solution 
returned  by  the  set  cover  based  algorithm  can  be  1.7  times  that  of  the  solution  returned  by  the 
Steiner  tree  based  algorithm.  These  two  algorithms  thus  offer  a  tradeoff  between  performance 
and  complexity. 


IV.  Simulation  Results 

In  this  section,  we  present  simulation  results  on  the  performance  of  the  two  approximation 
algorithms  (SCA  and  STA)  for  the  broadcast  problem  in  an  MR-MC  network.  We  consider  a 
dense  MR-MC  network,  where  all  the  nodes  are  within  each  other’s  transmission  range,  and  we 
aim  to  minimize  the  total  energy  consumption  of  the  broadcast. 

There  are  12  non-overlapping  channels  _/)  (1  <  i  <  12),  possibly  with  different  communication 
rates  rt,  available  for  the  MR-MC  network,  and  each  node  is  equipped  with  4  radios.  At  the 
beginning  of  the  broadcast,  each  node  randomly  selects  4  of  the  12  channels  for  its  4  radios. 
As  discussed  in  Sec.  I-C,  the  nodes  which  share  at  least  one  common  channel  form  a  clique, 
and  there  is  a  one-to-one  correspondence  between  the  cliques  and  the  faces  of  the  derived 
WSC.  The  weight  of  the  face  is  defined  as  the  energy  consumption  of  the  broadcast  within  the 
corresponding  clique,  i.e.,  the  sum  of  the  transmission  energy  and  the  reception  energy.  Let  S  be 
a  face  containing  k  +  1  nodes  and  {fsj  :  j  —  1,  2, ...,  q\  the  q  (1  <  q  <  12)  common  channels 
shared  by  the  k  +  1  nodes.  Assume  that  if  a  node  in  the  clique  is  selected  as  relay,  it  will  choose 
the  common  channel  with  the  maximum  communication  rate  to  transmit.  Then  the  weight  w(S) 
of  the  face  S  is  given  by 

w(S)  =  (. Ptx  +  kPrx) - - 

max 
.7=1, 

where  Ptx  and  Prx  are  the  transmission  power  and  the  reception  power,  respectively,  and  L  is  a 
constant. 

nA  random  simplicial  complex  A (n,D,p)  with  n  vertices,  dimension  at  most  D,  and  a  D-dimensional  probability  vector 
p  =  {pi,p2,  ...,pd}  is  constructed  in  a  bottom-up  manner:  first  n  vertices  are  fixed,  which  are  the  0-simplices  of  A,  and  then 
higher-dimensional  simplices  are  generated  inductively.  Specifically,  for  each  1  <  i  <  D,  after  all  the  simplices  with  dimension 
lower  than  i  have  been  generated,  consider  every  i- tuple  of  vertices:  if  they  have  formed  all  the  lower  dimensional  simplices, 
then  an  i-simplex  consisting  of  them  is  generated  with  probability  pi.  Notice  that  a  random  simplicial  complex  A(n,  l,p)  is 
the  random  graph  introduced  by  Erdos  and  Renyi  [37], 
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Fig.  7.  Average  total  energy  vs.  number  of  nodes.  Parameters:  Ptx= 1,  Prx  =  0.01,  L  =  100,  ry  =  i  for  1  <  i  <  12. 


In  Fig.  7,  the  average  total  energy  of  the  solutions  returned  by  SCA  and  STA  is  compared  with 
that  of  the  MST  with  respect  to  the  underlying  graph  of  the  WSC.  The  average  is  taken  over 
10  random  channel  assignments.  Notice  that  although  two  different  links  on  the  same  channel 
are  treated  separately  when  the  MST  is  derived,  the  transmission  energy  corresponding  to  them 
is  counted  only  once  to  exploit  the  wireless  broadcast  advantage  when  the  total  energy  of  the 
MST  is  computed.  We  see  that  the  performances  of  SCA  and  STA  are  extremely  close,  and  their 
performances  are  significantly  better  than  that  of  MST. 

V.  Conclusion  and  Future  Work 

In  this  technical  report,  we  study  the  minimum  cost  broadcast  problem  in  multi-radio  multi¬ 
channel  ad  hoc  networks,  where  the  total  cost  is  the  sum  of  the  costs  associated  with  the 
transmissions  during  the  broadcast.  We  formulate  it  into  a  fundamental  problem,  the  minimum 
spanning  problem,  in  simplicial  complexes.  Due  to  the  existence  of  the  simplices  with  dimension 
higher  than  edges,  this  minimum  spanning  problem  is  more  complex  than  its  counterpart  in  a 
‘conventional’  graph  for  which  simple  polynomial  algorithms  exist.  Specifically,  it  is  shown  to 
be  NP-complete  via  a  reduction  from  the  set  cover  problem.  We  thus  propose  two  approximation 
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algorithms  for  this  minimum  spanning  problem:  one  is  to  transform  it  into  the  connected  set 
cover  problem;  the  other  is  to  transform  it  into  the  node-weighted  Steiner  tree  problem  and  then 
apply  the  corresponding  algorithm.  Despite  their  distinct  approaches,  the  performance  of  both 
approximation  algorithms  is  shown  to  be  order-optimal.  Furthermore,  we  show  that  there  is  a 
tradeoff  in  terms  of  performance  vs.  complexity:  although  the  one  based  on  connected  set  cover 
has  lower  time  complexity  than  the  one  based  on  the  Steiner  tree,  the  performance  of  the  former 
is  sensitive  to  the  weight  range. 

As  a  starting  point,  we  have  assumed  that  the  channel  assignment  scheme  is  designed  indepen¬ 
dent  of  the  broadcast  scheme.  The  joint  optimization  of  the  two  schemes  will  further  reduce  the 
broadcast  cost.  Another  future  direction  is  to  develop  distributed  versions  of  the  approximation 
algorithms  for  the  minimum  cost  broadcast  problem. 
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